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Abstract 

We study a Fokker-Planck equation with double-well potential that is nonlocally driven by 

a dynamical constraint and involves two small parameters. Relying on formal asymptotics we 
identify several parameter regimes and derive reduced dynamical models for different types of 
phase transitions. 
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1 Introduction 



In this paper we investigate the different dynamical regimes in a Fokker-Planck equation with 
multiple scales that was introduced in [3] to describe the charging and discharging of lithium-ion 
batteries, a process that exhibits pronounced hysteretic effects [6, 2]. The model, to which we 
refer as (FP), governs the evolution of a statistical ensemble of identical particles and is given by 
the nonlocal Fokker-Plank equation 

TdtQ{x, t) = d:, (u'^d^Qix, t) + {H' {x) - a{t)) q{x, t)) . (FPi) 

Here H is the free energy of a single particle with thermodynamic state x G R, the probability 
density q{-, t) describes the state of the whole system at time t, and a reflects that the system 
is subjected to some external forcing. Moreover, r > is the typical relaxation time of a single 
particle and i/ > accounts for entropic effects (stochastic fluctuations). 

The model (FP) has two crucial features which cause highly nontrivial dynamics. First, the 
free energy H is a double-well potential, hence there exist two different stable equilibria for each 
particle. Second, the system is not driven directly but via a time-dependent control parameter. In 
our case this parameter is the first moment of q, that means we impose the dynamical constraint 
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xq{x, t) dx = i{t) , (FP2) 

where £ is some given function in time, and a direct calculation shows that (FP2) is equivalent to 

a{t) = [ H'{x)g{x, t) dx + Tl{t), (FP^) 
Jk. 

provided that the initial data satisfy xq{x, 0) dx = ^(0). The closure relation (FP2) imphes that 
(FPi) is a nonlocal and nonlinear PDE. Well-posedness was proven in [5] on bounded domains, 
but we are not aware of any result about the qualitative properties of solutions. 

An intriguing property of (FP) is that its dynamics involves three different time scales. On 
the one hand, there are the relaxation time r and the time scale of the dynamical constraint. On 
the other hand there is the scale of probabilistic transitions between different local minima of the 
effective energy Ha{x) = H{x) — ax, that means particles can move between the different wells 
due to stochastic fluctuations (large deviations) . Kramers studied such transitions in the context 
of chemical reactions [10] and derived the characteristic time scale 

^exp(^^^), (1) 

where Ai^o- is the minimal difference of energy between the local maximum any of the local 
minima of H„. In what follows we always assume that £ is of order 1, whereas both r and v are 
supposed to be small. 

Our goal in this paper is to identify different parameter regimes and to describe the asymp- 
totics of (FP) in the limit u,t ^ To this end we focus on strictly increasing constraints and 
describe four different mechanisms of mass transfer between two stable regions. The correspond- 
ing four types of phase transitions arc, roughly speaking, related to two main regimes, which we 
refer to as fast reacfAon regim,e and slow reaction region, respectively. The dominant effect in 
the fast reaction regime is mass exchange according to Kramers' formula. This appears for very 
small r and covers, as limiting case, also the quasi-stationary regime r = 0. The slow reaction 
regime, however, corresponds to very small v and Kramers' formula is not relevant anymore. In- 
stead, phase transitions are dominated by transport along characteristics and this causes rather 
complicated dynamics since localized peaks of mass can enter the spinodal region of H. 
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In both the slow reaction and the fast reaction regimes we are able to characterize the small 
parameter dynamics in terms of a few averaged quantities only. Detailed descriptions of the 
corresponding limit models are given in the introductions to Sections 2 and 3, respectively. 

1.1 Preliminaries about Fokker-Planck equations 

Before wc give a more detailed overview on the different dynamical regimes we specify our as- 
sumptions on H and review some basic facts about Fokker-Planck equations. 

1.1.1 Assumptions on the potential 

In this paper we assume that H is an even double-well potential that satisfies the following 
conditions, see Figure 1, 

x = X+ I 
^ /^-^ 



Figure \: H \s supposed to be a 'typical' but even double well potential, as for instance (2). The functions 
X_, Xq, and X+ denote the three branches of the inverse of H' . 

(Al) H is even, sufficiently smooth (at least C^), and H'{x) grows at least linearly as x ^ oo. 
(A2) There exist constants < < x** and < cr* such that 

(a) H'{x*) = —a* and H'{x**) = cr*, 

(b) -(7* < H'{x) < and H"{x) < for < x < x*, 

(c) -(7* < H'{x) < oo and H"{x) > for a; > x*. 

In particular, the inverse of H' has the three strictly monotone branches 

X_ : (—DO, a*] — > (—CO, — x*] , Xq : [—a*, a*] — )• [— x*, x*] , X_ : [—a*, oo) — )► [x*, oo) . 

(A3) The functions X+ o H' and X_ o H' are concave on the spinodal interval (— x*, x*). 

The assTimptions in (Al) and (A2) are made for convenience and might be weakened for the price 
of further technical and notational efforts; (A3) is a geometric condition that becomes important 
in the slow reaction regime that is discussed in Section 3. Notice that all assumptions are in 
particular satisfied for the standard double-well potential 

H{x) = (x2 - 1)^ (2) 

In what follows wc refer to (— oo, — x*) and (.x*, oo) as the stable intervals, whereas the spin- 
odal region (— x*, x*) is called the unstable interval. This nomenclature is motivated by the 
different properties of the transport term in (FPi). In both stable intervals adjacent characteris- 
tics approach each other exponentially fast, hence there is a strong tendency to concentrate mass 
into narrow peaks. In the unstable region, however, the separation of adjacent characteristics 
delocalizes any peak with positive width. 
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1.1.2 Thermodynamical aspects 



Fokker-Plank equations like (FP) are derived in [3] from First Principles and provide a thermody- 
namically consistent model for a many-particle system with dynamical constraint. In particular, 
the second law of thermodynamics can be stated as 

S = -V + ai, 

where 6 is the free energy of the system and T> the dissipation. They are given by 
V{q, a) = - [ -^({H'{x)-a)g{x) + u^d^g{x)Ydx>0 

'T JR Q\X) ^ ' 

and £ = % — where 

Hig) = / H{x)g{x)dx, S{g) = —v^ g{x) log g{x) dx 

JR JR 

denote the internal energy and entropy of the many-particle system, respectively. 

It is well known that the Fokker-Planck equation without constraint, that is (FPi) with a = 0, 
admits several interpretations as a gradient flow. There is for instance a linear structure, which 
has been exploited in [12] in order to derive the effective dynamics in the limit — t- 0. Of 
particular interest, however, is the nonlinear Wasserstein gradient flow structure, see [8, 9, 7, 1], 
as this structure is compatible with the constraint. More precisely, (FP) with ^ = is the 
Wasserstein gradient flow for £ on the constraint manifold g{x) dx = and f 7^ describes a 
drift transversal to this manifold. 

The entropic term 1/ is often supposed to be very small but it is important that u is positive. 
More precisely, without the diffusive term the qualitative properties of solutions would strongly 
depend on microscopic details of the initial data, and hence it would be impossible to characterize 
the limit r, — )■ in terms of macroscopic, i.e. averaged, quantities only, see [4]. A key observation 
is that the singular perturbation v'^d'^g regularizes the macroscopic evolution, at least for some 
classes of initial data, in the following sense. Microscopic small-scale effects are still relevant on 
the macroscopic scale, but they are independent of the initial details and affect the system in a 
well-defined manner. As a consequence we now obtain a well-posed limit model for macroscopic 
quantities. Another approach to ensure well-defined macroscopic behavior is investigated in [11]. 
The key idea there is to mimic entropic effects by assuming that each particle is affected by 
a slightly perturbed potential. The macroscopic evolution is then completely determined by 
the dynamical constraint, the macroscopic initial data, and the probability distribution of the 
perturbations. 



1.1.3 Dynamics in the unconstrained case 

We next summarize some facts about the dynamics of (FPi) with time- independent a and small 
parameters < z/, r ^ 1. For a ^ [— cr*, -|-cr*], the effective potential H^{x) = H{x)—ax possesses 
a single critical point x that corresponds to a global minimum. The system then relaxes very fast 
(on the time scale r) to its unique equilibrium state 




where Z is a normalization constant ensuring geqdx = 1. This equilibrium density ^eq has a 
single peak of width u located at x and decays exponentially as x ^ ±00. 
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The situation is different for a G (— f*, +0"*) since Hu now exhibits a double weU-structure 
with two wells (local minima) at x± = X±{a) that are separated by a barrier (local maximum) 
at xo = Xo{a). Initially the system relaxes very quickly, and approaches (for smooth initial data) 
a state composed of two narrow peaks located at the wells. Both peaks have masses m_ and m+ 
with m_ + m-i- = 1, but the precise values of m± depend strongly on the initial data. This fast 
transition reflects that each particle in the system is strongly attracted by the nearest well due 
to the gradient flow structure. 

The resulting state, however, is in general not an equilibrium but only a metastable state. 
The underlying physical argument is that particles can pass the energy barrier due to stochastic 
fluctuations. In the generic case, in which both wells have different energies, it is of course more 
likely for a particle to cross the barrier coming from the well with higher energy, and thus there 
is a net flux of mass towards the well with lower energy. This flux is, for small <C 1, given by 
Kramers' celebrated formula and guarantees that the system approaches its equilibrium on the 
slow time scale (1). The corresponding equilibrium solution is again given by (3) and has now 
two peaks with a definite mass distribution between the wells. Notice, however, that for small u 
almost all the mass of an equilibrium solution is confined to the well with lower energy. 

1.2 Overview on different types of phase transitions 

Due to the different time scales, the dynamics of (FP) can be very complicated, and we are 
far from being able to characterize the small parameter dynamics for all types of initial data 
and all reasonable dynamical constraints. We thus restrict most of our considerations to strictly 
increasing dynamical constraints and well-prepared initial; only in Section 3.3.2 we allow for 
non- monotone constraints. 

1.2.1 Monotone constraints and well-prepared initial data 

In what follows we consider functions i with 

e{0) < -x** , < c < i{t) <C <oo for all t> 0, (4) 

where c, and C are given constants. Since r is small, the system then relaxes very quickly to a 
local equilibrium state with a{0) = H'(i{0)) < —a*. We can therefore assume that the initial 
mass is concentrated in a narrow peak, that means 

q{x, 0) ^ 5i(^o){x) , (5) 

where the right hand side abbreviates the Dirac distribution at i{0)- An even better approxima- 
tion, that also accounts for the small entropic effects caused by < <C 1, is 

q{x, 0) ^ ^y^exp , « = H"m) > 0. (6) 

The dynamical constraint implies (7 > 0, so that the peak starts moving to the right and the 
system quickly relaxes to a new local equilibrium state. For sufficiently small times, that means 
as long as £{t) < — x*^,, the system can therefore be described by the single peak model 

g{x,t)^6e^t)ix), m_(t) = 1, m+(t) = 0, a{t) = H' {i{t)) , (7) 

where we write m_ (t) = 1 and m+ (t) = to indicate that all mass is confined in the left stable 
interval (— oo, — x*). Moreover, at some time > we have i{t) = x**, and for t > the system 
again relaxes quickly to a local equilibrium state, that means we have 

Qix,t)^5eit)ix), m_(t) = 0, m+(i) = 1, a{t) = H'im) , (8) 
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where m-{t) = and = 1 now reflect that all mass has been transferred to the second 

stable region (x*, oo). 

The key question is what happens between t = and t = t*. It was aheady observed in [3] that, 
depending on the relation between r and there are at least four types of phase transitions driven 
by rather different mechanisms of mass transfer from the left stable region into the right one. 
The main objective of this paper is to investigate the different regimes and to derive asymptotic 
formulas for the dynamics. 

1.2.2 Different regimes in numerical simulations 




Figure 2: Numerical solutions to the initial value problem (9) for several values of r and v. The curves 
Fstatc and Fphase from (10) are drawn in Black and Dark Gray, respectively; the Light Gray curve represents 
the graph of H' . 



We illustrate the different types of phase transitions for the driven initial value problem (4) 
and (5) by numerical simulations with 

H'{x) = X - 2 arctan {x) , l{t) = t-4, t e {0, 8) , (9) 

and initial data g{-, 0) as in (6). Figure 2 visualizes the numerical solutions by means of two 
curves 

rotate : t ^ Hit), y{t)), Tphase : t ^ {i{t), a,{m+{t) - m_(t))). (10) 

These curves represent the macroscopic state of the system and the (rescaled) phase field, respec- 
tively, and are defined by 

r r+oo pO 

y{t) = / H'{x)q{x, t) dx, m+{t) = / g{x,t)dx, m-{t) = / g{x,t)dx. 

Jr Jo J -oo 

The microscopic state of the system is illustrated in Figure 3 by snapshots of g at different times. 
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Figure 3: Snapshots of q at several times corresponding to four simulations from Figure 2; the shaded 

area represents the unstable interval. First row. A stable peak can pass through the unstable region. 
Second row. An unstable peak can split into two stable peaks. Third row. Two stable peaks exchange 
mass by a Kramers- type phase transition (cr const > 0). Fourth row. Two stable peaks exchange mass 
according to the quasi-stationary dynamics (cr w 0). 

We emphasize that Figures 2 and 3 also illustrate the different dynamical regimes for fixed 
< <C 1, in the sense that the limit r — >■ can be regarded as a passage from (^4) to (I). 
However, as our results will show, there is an exponential scale separation between the different 
regimes and thus it is very hard to capture all types of phase transitions in numerical simulations 
with the same value of v. 

The numerical simulations illustrate that there exist the following types of phase transitions. 

Type /, Example [A): At some time the narrow peak enters the unstable region due 
to the dynamical constraint and starts to widen due to the separation of 
characteristics. However, the transport is much faster than the widening, 
so that the peak can pass trough the unstable region. 

Type II, Examples {B), (C), (D), {E): The peak still enters the unstable region but 
now the widening is much faster than before. In particular, the unstable 
peak dclocalizes, and the system quickly forms new peaks in each of the 
stable regions. At a later time the left peak enters the unstable region, 
and the competition between transport and widening starts once more. 

Type III, Examples (F), (G): The peak does not enter the unstable region anymore. 
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Instead, at some position in the left stable interval the peak stops moving 
and starts loosing mass to feed another peak in the right stable interval. 
This Kramers-type process happens with a const > and goes on 
until the left peak has disappeared. 

Type IV, Example (I): This is the quasi-stationary limit. The phase transition is 
similar to Type III but happens with a 0. 

Notice that the types II and III imply hysteresis. In fact, if we revert the situation by driving 

the system with £ < and 1(0) > x**, the symmetry of the problem implies that the macroscopic 
state and the phase field are confined to the images of Tstate and Tphase under point reflections at 
0. 

1.2.3 Main results and organization of the paper 

Our main results are formulas for the macroscopic evolution in different parameter regimes. The 
corresponding scaling relations between r and u are summarized in Table 1 and the limit models 
are presented in the introductions to Sections 2 and 3. 

condition parameter regime type 

^ 1 — slow reactions I single-peak limit 



log 1/v 

a 

log l/u 



— i,,.*^ / < a < Ocrit slow reactions II piecewise continuous 

two-peaks evolution 

< p < I slow reactions OPEN PROBLEM 



T = |<p<oo fast reactions III limiting case of 

Kramers ' formula 

r = exp ^ 2^ < b < hcrit fast reactions III Kramers' formula 

T < exp ^ fast reactions IV quasi-stationary limit 

Table 1: Overview on the different scaling regimes for < t,v ^ 1. The constant /icrit is completely 
determined by the properties of whereas the constant Ocrit depends on both H and I. 

The rest of the paper is organized as follows. In Section 2 we derive asymptotic formulas to 
describe Type-/// transitions in the regime exp (— /icrit/^^^) *C r <C v^l"^ , where /icrit denotes the 
energy barrier of H. We first recall in Section 2.1 the formal derivation of Kramers' formula for 
the mass flux between the different wells of //o-- In Section 2.2 we then identify the critical value 
for a = a(T, u) for which such a mass flux can also take place in the constrained setting and 
show in Section 2.3 that small variations in a are sufficient to adjust the mass flux according to 
the dynamical constraint. Finally, in Section 2.4 we discuss Type-/y transitions as these can be 
regarded as limits of Kramers type transitions. 

In Section 3 we consider the scaling regime v = exp (—a/r) with a > and discuss Type-// 
transitions which contain, as a special case, also Type-/ transitions. In Section 3.1 we first neglect 
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all entropic effects and introduce a simplified two-peaks model that allows to understand how two 
Dirac peaks interact due to the dynamical constraint. It turns out that the dynamical constraint 
can stabilize a peak in the unstable region, but also that at some point a bifurcation forces 
both peaks to merge instantaneously. In Section 3.2 we introduce another simplified model that 
accounts for the stochastic fluctuations in the unstable region and allows to understand how an 
unstable peak delocalizes and splits into two stable peaks. In particular, we derive an asymptotic 
formula for the time at which such a splitting event takes place and introduce the mass splitting 
problem that determines the mass distribution between the emerging peaks. In Section 3.3 we 
finally combine all result and characterize the limit dynamics as intervals of regular transport 
that are interrupted by several types of singular event. 

2 Fast reaction regime 

In this section we show that Kramers type phase transitions are also relevant in presence of the 
dynamical constraint as long as 



where /icrit = ^(0) — minajgjj H{x) denotes the energy barrier of H. The key idea is that a 
Kramers type phase transition occurs during a time interval (^1,^2) in which a is positive and 
almost constant. During this time interval a only changes to order 0(z^^) but this is sufficient to 
accommodate the dynamical constraint. Kramers' formula therefore allows to understand phase 
transitions of type III, and hence that a stable peak suddenly stops moving and starts loosing 
mass to feed a second stable peak. 

The situation is different for r <C exp (— /icrit/^^^) since then we expect to find phase transitions 
of type IV, that means the mass flows towards the second well as soon as it is energetically 
admissible. This regime is governed by the quasi-stationary approximation but can also be 
regarded as a limiting case of Kramers regime. 

Our main result concerning the fast reaction regime combines the formal asymptotics for the 
Kramers regime and the quasi-stationary approximation and can be stated as follows. 

Main result. Suppose that the dynamical constraint and the initial data satisfy (4) and (5), and 
that T and v are coupled by 






for some constant b G (0, /icrit)- Then there exists a constant ai, G (0, cr*) such that 



1. the dynamical multiplier satisfies 




where ti and t2 are uniquely determined by i{ti) = X_(a"(,) and £{t2) = X^{ab) 



2. the state of the system satisfies 
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where = 1 — m^{t) and 

' 1 

m-{t) = < 



X+{ah)-m 







for t <ti, 
for ti < t < t2 , 
for t>t2. 



Moreover, the assertions remain true 

1. with (76 = ifT< exp (^-^) , 

2. with ab = (T* if T <^ but r > exp (— ^) for all b > 0. 

To justify the limit dynamics we review Kramers' argument for constant a in Section 2.1. In 
Section 2.2 wc then derive similar asymptotic formulas for the constrained case, which allow us to 
adjust the mass flux according to the dynamical constraints in Section 2.3. Moreover, in Section 
2.4 we discuss the quasi-steady approximation, which governs the regime < r ^ exp (— /icrit/z^^) • 

We finally mention that the limit energy is given by 

E :=m-H{X^{a)) +m+H{X+{a)), 

and evolves according to 



E = al-DbX{a=a,}i Db-.-- 



H,,{X_{a,)) - H^,{X+{a,)) 



>0, 



X+{ab) - X_{ab) 

where H^r is defined by H^{x) := H{x) — ax and X{(T=ab} denotes the usual characteristic function. 
2.1 Kramers' formula in the unconstrained case 

To derive Kramers' formula for the unconstrained case we consider the Fokker-Planck equation 
(TPi) with fixed a G (— o"*, cr*) and use the abbreviations 

Xi:=Xi{a), ai := \H" {xi)\ , ie{-,0,+}. 

At first we approximate g outside the local maximum xq = Xo(cr) for small u by the ansatz 



kHa{x) = H(x) - ax 




Figure 4: Particles can cross the energy barrier between the different wells of due to random fluctua- 
tions, and this gives rise to an effective mass transfer from the well with higher energy towards to the well 
with lower energy [Kramers type phase transition). 
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q{x, t) 



( HAx)\ „ ^ 
/x_(rjexp I 2 — ) lor x < xq, 

exp I 2 — I tor X > xq. 



This is the outer expansion and reflects the assumption that the system has a peak in either of 
the stable regions, where the masses are given by 



m±{t) := ±//±(i) j exp ^ — j '^^ 



For small u we can simplify the integrals using Laplace's method, that means we expand H^r 
around x± to find 

m±[t) « — A*±(ij exp I — I . 



The mass exchange between both peaks is then determined by 

rm±(t) = Ti?(t), (13) 

where R{t) is the mass flux at .xq. The key idea behind Kramers' formula is that can be 

computed from the quasi-stationary approximation of q near xq. More precisely, with the change 
of variables y = (x — xq) / v we approximate 

v'^dxQit, Xq + vy) + H'„{xq + i^y)g{t, xq + vy) ~ ~R{t), 

and obtain the inner expansion 

Q{t, xo + vy) « exp 5 C{t) / exp ^ h 



where C{i) is a constant of integration. For small v and y S> 1 we can simplify the integrals by 
Laplace's method to obtain 



/ 



\^ (^M?^) dS . fexp di = /Xe.p 



'0 V / ^0 V 

and using a similar formula for y <C — 1 we flnd 



Q{t, xo±v \y\) « exp -2 J [C{t) T — exp \^—^ j j ■ (14) 

In order to match the outer and the inner expansions, we consider x ^ xq and compare the 
asymptotic formulas (12) and (14). Since both contain the factor exp [—Hc^{x)/v^^, we equate 
the time dependent coefficients and arrive at the matching conditions 

^^(t) = C(t)T«,/;^exp^^^(^°) 



V Y 2ao ^ V 

We finally eliminate C{t) and find the desired expression for Kramers' mass flux, namely 



'2tt 

271 



Rit) ^ ^(^-(t) - ,.{t)) exp (-^) (15) 

m_(t)-y/Q;_aoexp I — — I — m+(t)-y/a+Q!o exp I — — I I . 



1 / H^{xo) 
exp 5 — 
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Combining this with (13) we easily verify that the characteristic time for Kramers mass transfer 
is given by 

/ Ha{xo) -Taax{H^{x+),H^{x-))\ 
rexp 5 . (16) 



We also notice that in the generic case H„{x+) ^ Hu(x^) and for small z/ the mass transfer is 
essentially unidirectional on the time scale (16), that means the mass flows from the well with 
higher energy to the well with smaller energy, see Figure 4. 

2.2 Kramers' formula in the constrained case 

We now derive a self-consistent description for Kramers type phase transitions in the presence of 
the dynamical constraint. To this end it is convenient to replace r by the parameter b defined in 
(11), and to consider the functions 

h±{a) := H^{Xo{a)) - H^{X±{a)) = H{Xo{a)) - H{X±{a)) + a{X±{a) - Xo{a)) . 

These functions are well-defined for \a\ < a* and satisfy dh-/ da < < dh+f da with 

h-{—a^) = /i+(cT*) = H{x^) - H{x^^) + a*{x^ + x**) =: h^, h+{—a^) = h-{a^) = , 
and /i_(0) = /i+(0) = -ff(O) — min-rgK -ff(x) =: /icrit < 

In terms of 6 and h±, the rescaled flux R{t)/T from (15) can be stated as 

^ ^ m_(t)r_(a(i)) - m+{t)r^{a{t)) , (17) 

where 



r^{a) := ^"^^""^"^ exp ( ^-^^^^ ] . (18) 



and ai{a) := \H" {Xi(a))\ as above. 

We next present some heuristic arguments for the dynamics of the rescaled flux terms m_(t)r_(t) 
and m+{t)r+{t). To this end we assume 

< 6 < /icrit, 

and recall that, due to our assumptions on the initial data, the system evolves for small times 
according to the single peak evolution (7). In particular, the peak reaches the critical position 
— .T*H< at time to with o-(to) = — o"*, and the dynamical constraint implies a{to) > 0. Assuming 
that a changes regularly at to, we then conclude that 

h-ia{to)) = K > 0, -^h^iait)) <0, h+{aito)) = 0, -^^h+{a{t)) > 0. 

Consequently, for small u and sufficiently small times t > to we expect to find r_(t) <C 1 <C ^^(t), 
so crossing the energy barrier is very likely for a particle in the right well but very unlikely for 
a particle in the left well. However, the net transfer across the energy barrier is very small since 
there are essentially no particles in the right well. We thus expect that the partial masses m± stay 
constant in the limit — t- 0, so the system can still be described by the single peak approximation 
(7). At some later time to > to we have a (to) = and the fluxes r_(to) and r+(io) have the same 
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order of magnitude. However, both are very small due to b < /idit = h± (ioj , and so there is, for 
small v, still no effective mass transfer between the two wells of H^^. 

The situation changes completely at time ti defined by a{h-{ti)) = b. At this time, r_(ti) 
becomes sTiddcnly of order one and wc can no longer neglect particles that move from the left 
well to the right one. The other flux r+, however, is now very small as it is very unlikely that a 
particle moves the other way around. 

As explained above, the main idea in the dynamical case is that there exist a constant ab and 
a time t2 such that a{t) ~ (jfy for all ti < t < t2- Kramers' mass flux can hence stabilize to 
continuously transfer mass from the left well to the right one. At time t2, all mass has been 
transferred to the right well and the system again evolves according to the single peak evolution, 
now given by (8). 

Before we describe the details of the mass transfer we proceed with two remarks. First, the 
above assumption b < /icrit is truly necessary: For /lent < b < h^^ there is still a time ti with 
h-(a{ti)) = b, but then we have h+{a{ti)) < h-{a{ti)) and hence r+(ti) » r_(ii), which shows 
that a net transfer from the left well to the right one is impossible. Moreover, for b > h* both r_ 
and r-)_ are always very large. In both cases we expect that the phase transition occurs when o" ~ 
and is not governed by Kramers' formula anymore but by the quasi-stationary approximation. 

Second, the mass flux is already determined by the dynamical constraint and the assumption 
a{t) ~ (Jfo. In fact, the constraint (FP2) implies that 

i{t) ^ X_{ab)m-{t) + X+{ai,)m+{t) 

and thus, since a does not change much, 

£{t) « X-{ab)m-{t) + X+{ab)m+{t) = - X+{ab))m-{t) . 

As a consequence we obtain 

for t2 < t < ts, where ts > t2 is defined by m_(t3) = 0, and using £{t2) ~ X-{ab) it is easy to 
check that £{ts) ^ X+{ab)- 

2.3 Adjusting the mass flux by small variations of the multiplier 

To identify the formulas that relate the mass flux self-consistently to small temporal changes in a 
we consider only times t with ti < t < t2 and assume that both m_ and m+ are strictly positive. 
Of course, in order to match the resulting approximations for g to the single-peak evolution for 
t < ti and t > t2 we must introduce transition layers at t ~ ti and t ^ t2 corresponding to 
m_ « 1 and m_ « 0, respectively, but since these transition layers do not contribute to the limit 
model, we do not investigate them in detail. 

Case 1 : limj^^o ^ > 0. The critical value cTb is defined by 6 = h-{ab)- Thanks to 6 < /icrit we 
find < (Tft < (J* and h-{ab) > /icrit > h^{ab), so (18) yields 

^ ( h.{ab)-h_{aim ^^^^^ ^ ^ ^2^^ 

ZTT \ J 

Since r_ (t) must be of order one, we introduce the rescaled multiplier 
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and simplify (20) by expanding h-{a) around ab- Using (13) and (17) we then conclude that 
Kramers' formula implies the mass transfer law 



-m_(t) - exp (|^'_(a,)| ^(t)) 



Comparing this with (19) we finally conclude that Kramers type phase transitions comply with 
the dynamical constraint if and only if 



In 



27r i{t) 



^^'^ " KM ■ ^ ^ 

This is the heart of our argument. If ijj evolves according to (21), then a — o";, is of order i/^, and 
Kramers' formula provides a mass flux that satisfies the dynamical constraint. 

Case 2 : lim^_»o 6 = . In this limiting case, the phase transition happens when a is close to cr*, 
so both a;_ and xq are close to —x^. Moreover, the constants a_ and approach zero, and 
hence we can no longer use the asymptotic expressions from the first case. However, if we expand 
all relevant quantities around cr* and — x*, it is still possible to derive an asymptotic formula for 
Kramers' flux that is consistent with the dynamical constraint. 

Thanks to the identities X_((T*) = Xo(o-*) = — x*, H"{—x*) = and 7 := —H"'{—x*) > 0, 
we deduce from the definition of X_ and Xq that 

X_(a) = -X* - + 0{a. - a), Xo{a) = -x, + ^^^^^^ + 0{a, - a) . 

To leading order in cr* — cr, we therefore find 

^/ a-{a)a+{a) « -^27(0-* - a) , 

as well as 

H^{Xo{a)) - H^{X-{a)) ^ -^(^{Xo{a) + x^f - {X^{a) + x^f) + (a* - a){Xo{a) - X+{a)) 

= ^(a* - a){Xo{a) - X+(a)) « ^(a* - af, 



so (15) can be simplified to 



o/.x rn-{t)^ 1/2 / 4^/2 [a^ - a{t)) 

^(*) ~ ./^ - exp ( - — ^ ^2 1 • (22) 



3/2' 



In Kramers' regime this flux R{t) should be of order r. On the other hand, the asymptotic formula 
(15) holds only if Hfj{xQ) — H„{x-) 2> v^, and thus we shall guarantee that v~'^{a^ — cr)^/^ ^ 1. 
Both conditions can be satisfled if 

In fact, if we define for given v ^1 and r ^ v"^/^ the large parameter K by 

TT \ 3^ / 
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then R{t) becomes of order r if (T* — cr is of order K'^v*/^ ~ i^v"^ In {l/v)^^^^ . 

Finally, we proceed as in the first case. Inserting (22) into (13) gives a formula for 
m_(t)/m_(t) in terms of cr* — a{t). On the other hand, the dynamical constraint implies 

m_(t)/m_(t) ^ -— / I 



'2iX^ j y 

via (19), and eliminating m_(t)/m_(t) we obtain cr* - (T(i) = 0(^^2/4/3) jn terms of ^(t). 
2.4 Phase transitions in the quasi-stationary hmit 

To conclude this section we show that the quasi-stationary approximation of (FP) describes 
phase transitions with a{t) « 0. Notice that such Type-IV transitions can be regarded as limits 
of Type-Ill transitions in the sense that cxft \ as 5 /icrit- 

In the quasi-stationary limit we approximate g by the equilibrium solution (3) that corresponds 
to the current value of cr. In other words, we set 



exp 

g{x, t) ^ ea{t){x), Qaix) := 



H{x) + ax 



i/2 



Zia) 

where the normalization factor Z{a) := Z^{a) + Z+{a) is given by 

^ / X f^"^''^ ( il{x)^ax\ , 

Z±{a) := ^ exp [ J dx . (23) 

The dynamical multiplier (T(t) is then determined by the dynamical constraint via 

^(i) = / Qa{t){x)^x. 

JR 

We now derive asymptotic formulas that characterize the quasi-stationary dynamics for small v. 
At first we notice that Laplace's method applied to (23) with a{t) ^ yields 

Z±{(T) « , , . exp 2 • 

V«±(^) V / 

and hence 

^ ^ I Sx+{a){x) for (T<0. 

In particular, the system evolves according to the single-peak approximation as long as a has a 
sign, and a phase transition can occur only for a = 0. To describe the details of such a transition 
it is convenient to rescale a by a = u'^ip. 
Using the expansion 

H{X±{a) + vy) - a{X±{u) + vy) = H{X±{a)) + lH"{X±{0))iyY - ^X±{0) + 0{a\au,u^) 
and employing Laplace's method once more, we find 



V«±(o) V / 
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and hence 



m_(t)/m+(t) ~ exp ( - 2X+(0)V'(i)) ■ 



On the other hand, the dynamical constraint, see (19), provides 



m-{t) Ri 



X+{Q)-i{t) 
2X+(0) 



2X+(0) 



and we conclude that the rescaled multiplier evolves according to 



m = - 



In 



) 



2X+{0) 



Notice that this formula is well defined for all times t with ti < t < t2, where ti and t2 are 
defined by i{ti) = X_{0) = -X+{0) and e{t2) = +X+{0), and satisfy limt^ti^{t) = -oo and 
limj^tj = +00 • 

3 Slow reaction regime 

This section concerns the effective dynamics of (FP) in the slow reaction regime: Both r and 
v are still supposed to be small, but v is so small that 'reactions', that means continuous mass 
transfer between the stable regions as described by Kramers' formula, are not relevant anymore. 
Instead, the dominant effect in (FP) is now transport along characteristics and this gives rise 
to new phenomena. In particular, localized peaks can enter the unstable region and peaks can 
split or merge rapidly. Notice, however, that the small entropic effects caused by > are still 
relevant and cannot be neglected. They guarantee that each peak entering the unstable region is 
basically a rescaled Gaussian, and hence that such a peak behaves in a well-defined manner. 

We now introduce an informal concepts that is motivated by numerical simulations and turns 
out to be useful for describing the asymptotic dynamics in the slow reaction regime. We say the 
system is in a two-peaks configuration if there exist two positions xi,X2 and two masses mi,m2 
with < 1712 = 1 — mi < 1 such that the state can be approximated by 



A two-peaks configuration is called stable-stable if xi < —x^ and X2 > x^.., but unstable- stable if 
—x^ < xi < x^ and X2 > x^. Moreover, in case that one of the masses vanishes, we refer to a 
two-peaks configuration as a {stable or unstable) single-peak configuration. 

Our main result in this section is an asymptotic description of the dynamics in a certain param- 
eter regime for r and v. The corresponding limit model is illustrated in Figure 5 and can be 
summarized as follows. 

Main result. Suppose that the dynamical constraint and the initial data satisfy (4) and (5), and 
that T and v are coupled by 



for some constant a > 0. Then, in the limit r — >■ the dynamics can be described in terms of 
single-peak and two-peaks configurations according to one of the following scenarios. 



g{x, t) ^ mi(t)(5^^(t)(x) + m2{t)6^^(^t){x) . 




(24) 
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stable transport (in x < —a;*) 
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unstable transport 
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Figure 5: Flowchart for the hmit dynamics in the slow- reaction regime with strictly increasing dynamical 

constraint. Intervals of quasi-stationary transport (of cither a single-peak or a two-peaks configuration) 
are interrupted by several types of singular times (corresponding to switching, merging and splitting of 
peaks). Gray and Black arrows indicate Type-I and Type-II phase transitions, respectively. 



1. Transport with splitting (Type-II phase transition).- There are intervals of quasi- stationary 

transport which are interrupted by singular times. More precisely, depending on i and the 
value of the parameter a there exists an integer K > 1 such that there are K + 1 switching 
times, K splitting times, and a final merging time. 

(a) During the intervals of quasi- stationary transport, the peaks do not exchange mass and 
are move according to the dynamical constraint. 

(h) At each switching time a stable peak reaches the position — to become unstable; each 
switching time is followed by a splitting or the merging time. 

(c) At each splitting time an unstable peaks splits and its mass is instantaneously trans- 
ferred to the stable regions. In particular, the system jumps from an unstable- stable 
two-peaks configuration (or, initially, from an unstable single-peak configuration) to an 
emerging stable-stable two-peaks configuration. The precise values for the masses and 
the positions of the emerging peaks are determined by a mass splitting problem. 

(d) At the final merging time the two peaks of an unstable- stable configuration merge - 
either continuously or dis continuously - to form a single stable peak located in the 
region x > x^. 

2. Pure transport (Type-I phase transition); The system is always in a single peak configura- 
tion and there exist only two singular times corresponding to switching ( the peaks enters the 
unstable interval) and trivial continuous merging (the peak leaves the unstable interval). 

Notice that it is practically impossible to perform numerical simulations with u exponentially 
small in r and r much smaller than 1. The phase transitions of type I and type II presented in 
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Figures 2 and 3 are therefore not directly covered by our limit model, but a close look to the 
numerical data reveals that they are likewise dominated by the interplay between transport and 
widening of unstable peaks. It remains a challenging task to derive next order corrections to 
replace the singular times by transition layers with width depending on r. Another interesting 
but open question is whether there exist scaling laws different from (24) that give rise to other 
reasonable slow reaction limits. 

To justify the limit dynamics we introduce several reduced models that allow us to study each of 
the different phenomena in a simplified setting. In Section 3.1 wc derive a two-peaks model that 
describes the transport of two peaks in terms of a simple ODE system. Moreover, this model also 
reveals that the separated peaks in an unstable-stable configuration can merge instantaneously 
due to the dynamical constraint. In Section 3.2 we investigate the splitting of unstable peaks. To 
this end we propose a peak-widening model that accounts for the small cntropic effects and allows 
us to derive a deterministic equation for the width of an unstable peak. In particular, we show 
that for small r it can happen that this width blows up almost instantaneously, which in turn 
gives rise to rapid mass transfer from the unstable towards the stable regions. Wc then discuss 
the mass splitting problem, which consists of solving a nonlocal transport equation in order to 
determine how much mass is transferred to each of the stable regions. Finally, in Section 3.3 
we combine all partial results and characterize the limit dynamics of the original model (FP). 
In particular, we derive explicit formulas for the iterative computation of all switching, splitting, 
and merging times. Moreover, at the end we sketch the slow-reaction limit for non-monotone 
dynamical constraints. 

3.1 Two-peEiks approximation: Transport and merging of peeiks 

A major tool in our analysis is a simple two-peaks model (TP) which describes the essential 
dynamics of (FP) as long as the system is a two-peaks configuration. The model governs the 
evolution of the peak positions xi and X2, and reads 

TXy2=<^-H\xi/2), (TPi) 
£ = miXi -\- 1712X2 (TP2) 
a = miH'{xi) + m2H'{x2) + t(. . (TP^) 

Here mi and m2 = 1 — mi denote the constant masses of the peaks. Notice that (TPi) is just 
the characteristic ODE for (FPi) with z/ = 0, and this implies that (TP) is also a constrained 
gradient flow corresponding to the energy 

E := miH{xi) -\- m2H{x2) ■ 

In particular, using the dissipation 

D := Tmix\ + Tm2x\ = — ^ (if'(xi) — a)^ H ^ {H' {X2) — cr)^, 

the energy balance is given hy E = —D -\- at. 

Since r is small it seems natural to neglect the time derivatives in (TPi) and (TP2). This gives 
rise to the quasi- stationary approximation to (TP), which consists of the algebraic equations 

H'{xi) = H'(x2) = <T, mixi -|- m2X2 = L (25) 

For our analysis it is important to understand in which sense (25) approximates (TP). This 
problem is not trivial because the non-invertibility of H' implies that (25) has multiple solutions 
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Figure 6: Left panel: Numerical solution to the two-peaks model (TP) with < r ^ 1. As long as 
t — to is sufficiently small, each peak is located in one of the stable regions, but at time ti the peak at xi 
enters the unstable region (switching configuration A). At a later time t2 the quasi-stationary two-peaks 
approximation ceases to exist, and the system jumps almost instantaneously from configuration B to the 
single-peak configuration C (discontinuous merging of peaks). Right panel: Cartoon of all steady states 
in the (xi, a:;2)-plane. The solution from the left panel is initially confined to the curve H'(x\) = H'{x2) 
with X2 > Xl but jumps to the diagonal Xi = X2 at time <2- 

for \i\ < X**, and thus we have to understand which solution branches are dynamically selected 
by solutions to (TP). 

In the limit r — )■ 0, solutions to (TP) exhibit two important dynamical phenomena which cor- 
respond to changing the solution branch of (25). Both phenomena seem to be counter-intuitive 
at a first glance but are a consequence of the dynamical constraint. They can be described as 
follows. 

First, the dynamical constraint can drive the system from a stable-stable configuration to an 
unstable-stable configuration, that means the stable peak at X2 can stabilize an unstable one at 
Xl due to the dynamical constraint. The second effect is that this stabilization can break down 
eventually. When this happens, the separated peaks merge almost instantaneously to form a 
single stable peak. 

To describe both phenomena we consider times t > to and suppose that the dynamical con- 
straint £{t) is smooth and strictly increasing. We also suppose that the initial data for (TP) are 
well prepared via 

xi{to) = X_{a{to)), X2{to)=X+{a{to)), e{to) = miXi{to) + m2X2{to) ■ (26) 

To elucidate the key ideas we proceed with discussing some numerical results and present some 
semi-rigorous analytical considerations afterwards. 

3.1.1 Dynamics of the two-peaks model 

The left panel of Figure 6 depicts a typical numerical solution to (TP) with < r ^ 1, potential 
(2), and initial data as in (26). The simulation reveals the existence of two critical times ti and 
t2 that separate the three different regimes 

Ri : to < * < —X** < Xl < — x*, 

R2 ■ tl < t < t2, —X* < Xl < -l-x*, 

R3 : t2 < t < ts, -l-x* < Xl < -l-x**, 
where t^ denotes the final simulation time. For all times < t < ^3 we have 

H\xi) Ri Cr, H'{x2) ~ (T, X* < X2 < X**, 
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and thus we expect that in the hmit r — )■ the system can in fact be described by quasi-stationary 
peaks. The details, however, are different for i?2, and i?3. More precisely, initially we have 

Ri : xi{t) ^ X_{a{t)), X2{t) ^ X+{a{t)), 

that means the solution resembles a stable-stable configuration, and i > implies xi > 0, X2 > 0, 
and & > 0. At time ti, which is defined by xi(ti) = — x*, the peak at xi enters the unstable 
region and the configuration switches to unstable-stable. This means 

i?2 : xi{t) ^ Xo{a{t)), X2^X+{a{t)), 

and hence xi > 0, X2 < 0, and < 0. 

At the second critical time t2, the two-peaks approximation breaks down, that means the 
system can no longer be described by two separated peaks. Instead, both peaks merge discon- 
tinuously, in the sense that the system jumps from an unstable-stable two-peaks configuration to 
a stable single-peak configuration. The evolution for t > t2 is still quasi-stationary, but involves 
only a single peak that evolves according to 

i?3 : xi{t) ^ X2{t) ^ lit). 
3.1.2 Failure of the quasi-stationary approximation 

As illustrated in the right panel of Figure 6, the existence of the critical time t2 can be understood 
as follows. The quasi-stationary two-peaks approximation imposes the constraints 

H'{xi) = H'{X2), Xl < X2, X^ <X2 < X**, 

which define a smooth curve in the (xi, X2)-plane. This curve has the two branches 

S_ := X+{u)) : \u\ < a,], Bo := {{Xo{a), X+{a)) : \a\ < a*}, 

which meet smoothly at the point (— x*, x**). Due to the dynamical constraint the system is 
further confined to 

:= |(xi, X2) : miXi 777,2X2 = ^1 , 

which is a straight line with slope —mi/ 1712- It can now easily be seen that in the quasi-stationary 
approximation for to < t < ti the state of the system corresponds to the unique intersection point 
of i?_ and G^, which moves towards (— x*, x**) since £ is increasing in time. At time ti, the system 
crosses this point (— x*, x**), and for sufficiently small times t > t\ the state of the system is 
given by the unique intersection point of Bo and G^. At time t2 however, the line becomes 
tangential to Bq and the system can no longer follow the curve Bq due to £ > 0. Instead both 
peaks merge rapidly and the system jumps almost instantaneously to {I, £), which is the only 
intersection point of G^ and the diagonal xi = X2. 
The tangency condition for the merging reads 

miH'\x2) + m2H"{xi)=0, 

where the left hand side is positive for all times to < t < t2- Of course, if 1712 is very small, then 
the slope of G^ is very negative and the tangency condition cannot be satisfied. In this case, 
there is a continuous transition from the unstable-stable two-peaks configuration to the stable 
single-peak configuration since both peaks merge via xi = X2 = i. It is then natural to define t2 
as the time at which this continuous merging takes place. 
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The above considerations also apply to the limiting case mi = 1 and m2 = 0, provided that we 
accept that X2 is undefined. In this case, the quasi-stationary approximation reads xi{t) = i{t) 
and the switching time ti is defined by i{ti) = — x*. At time t2 with i{t2) = +a;* the single peak 
leaves the unstable interval and in view of the above discussion it makes sense to interpret this 
event as (trivial) continuous merging. 

We emphasize that our characterization of the unstable-stable evolution relies on condition (A3), 
that means on the concavity of the function xi G [— .t*, ^ X2 = X^{H' (xi)). Otherwise it 
may happen at time t2 that the system does not jump to the diagonal but instead to another 
unstable-stable configuration. 

We also mention that discontinuous merging of peaks implies that a jumps down to H'{£), 
and that the emerging single peak is stable due to x^ < i < x^^. Moreover, the energy of 
the system also jumps down. This is in line with the gradient flow structure of (TP), and 
holds because the energy decreases along the straight line that connects (xi, X2) to (i, £). In 
fact, with ^(A) = miH{{l - X)xi + M) + m2H{{l - X)x2 + X£) we find dE/dX < thanks to 
£ = mixi +m2X2, xi < X2 and H'[{1 - X)xi + X£) < H'[{1 - X.)x2 + A^). 

3.1.3 Stability of unstable-stable two-peaks configurations 

We finally show that the quasi-stationary approximation is linearly stable until it ceases to exist. 
To this end, we consider a given solution to (TP) with < r ^ 1 and initial data as in (26), and 
denote by xi and X2 the quasi-stationary approximation to (TP). This reads 

xi{t) = X_/Q{a{t)), X2{t)=X+{a{t)), mixi{t) + m2X2{t) = £{t), (27) 

where the definition of xi involves X- for to < t < ii but Xq for ti < t < t2- Here, t2 is the 
time at which the quasi-stationary approximation ceases to exist, and ii is defined by ^(ti) = 
—mix* -|- m2X** and denotes the switching time at which the quasi-stationary approximation 
enters the branch Bq. 

Making the ansatz Xi{t) = Xi{t) + Tyi{t), we easily derive the linearized model 

ryi = m2(^H" {x2)y2 - H"{xi)y^ + 51, 
ry2 = mi [H"{xi)yi - H"{x2)y2^ + 92, 
a = a + T(miH"{xi)yi + m2H"{x2)y2 + ^ 

with 9i = £ — Xi- By construction, we have miyi + m2y2 = and hence 

Tyi = -Cyi + 9i, 

where C = m.\H"{x2) + m2H" {xi). Notice that C(^); 9i{t) and g2{t) depend continuously on 
t since the function £ is smooth. The Variations of Constants Formula now reveals that the 
quasi-stationary solution is dynamically stable as long as > 0, that means as long as t < t2- 
Moreover, a similar analysis reveals that the quasi-stationary single-peak solution xi = X2 = ^ is 
dynamically stable provided that H"(£) > 0, which is satisfied for t > i2 due to ^(^2) > a;*. 

3.2 Entropic effects: Widening and splitting of unstable peaiks 

In the previous section we have seen that the two-peaks model (TP) allows a stable peak to enter 
the unstable region. Our goal in this section is to show that under some conditions the same 
is true for the original model (FP), but also that the entropic terms trigger new phenomena. 
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To point out the main idea we start with some heuristic arguments. Afterwards we derive and 
investigate a simpUfied model that ahows to understand the key effects in a more formal way. 

Suppose that we are given a solution to (FP) that at time to consists of two narrow stable peaks 
located at positions —a;** < xi{tQ) < — and < X2(to) < Suppose also that the width 
of these peaks is sufficiently narrow such that |<T(to)| < c*- For small r we then expect that the 
system relaxes very fast to a meta-stable state of (FP) that meets the constraint £{to). Without 
loss of generality we may hence assume that the initial peaks have width of order v and that the 
initial data are well prepared in the sense of (26) 

For sufficiently small times t > to we can expect that (FP) follows the two-peaks model (TP), 
that means both peaks have width of order v and are located at stable positions xi G (— x**, — x*), 
X2 G (x*, X**) with xi > and ±2 > due to £ > 0. At some time ti, however, the peak located 
at xi reaches the critical position —x^,, and this time can be estimated by i'(ti) ~ — ?7?,ix* + m2a;**. 

When the first peak has crossed — x*, its width widens very quickly because the characteristics 
of the transport term now separate exponentially with local rate —H"{x)/t. However, since the 
width of the first peak is initially exponentially small in r, it remains small for some times 
although it is surely much larger than u. Moreover, the second peak located at X2 still has width 
of order i/ as it remains confined to the stable region x > x*. Combining both arguments we 
conclude that the system can be approximated by (TP) even at times t > ti provided that t — ti 
is sufficiently small. The condition ^ > then implies xi > 0, X2 < and a < 0. 




Figure 7: Schematic representation of the shape of g at different times (for better illustration with 
respect to the graph of H', that means the plots show g{x, t) + a{t) against x). Due to the dynamical 
constraint, the peak located at xi switches from stable (A) to mistablc (B) at time ti, and its width starts 
to grow exponentially. At time t2, the width of the peak becomes suddenly of order 1 (C) and the mass 
splitting problem describes that the system jumps almost instantaneously to a new stable-stable two-peaks 
configuration {D). 

The key question now is how long the first peak remains localized. Depending on the scaling 
parameter a and the mass distribution between the peaks, it can happen that the first peak 
remains localized until both peaks merge continuously via xi ~ 2:2 ~ ^ ~ x*. In this case the 
first peak can in fact pass through the whole unstable region. It can also happen that the peak 
remains localized till the quasi-stationary two-peaks approximation ceases to exist. In this case 
both peaks merge instantaneously and discontinuously to form a single stable peak, but (FP) still 
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behaves like (TP). 

There is, however, a third possible scenario, in which the width of the first peak becomes of 
order one before both peaks can merge continuously or discontinuously. Below we will argue that 
if this happens at all, it happens at a precise time t2- More precisely, we show that there is a 
time t2 at which the width of the first peak blows up instantaneously for small r. Some amount 
of the mass of this peak is then transported along characteristics to the left until it creates a 
new peak in the stable region x < — x*. The remaining part, however, is transported towards the 
other stable region x > to feed the second peak. Since the transport along characteristics if 
very fast, we expect that in the limit r — t- the first peak splits and disappears instantaneously 
and that the system jumps to another stable-stable configuration. After this jump the new peak 
in the stable region x < — a;* has mass smaller than mi, and this implies that the mass of the 
second peak is larger than 1712- 

3.2.1 A simplified model 

To analyze the relevant phenomena we study a simplified peak-widening m,odel (PW), which 
approximates the second peak at X2 by a Dirac mass but keeps the probabilistic description for 
the first peak located at x\. This gives rise to the equations 

rdtQ = d^(v''d:,Q + {H'{x) -a)Q^, (PWi) 

TX2 = Cr- H'{X2) , (PW2) 

I = m\ I xgdx + m2X2, (PW3) 
Jr 

a = mi H'{x)Qdx + m2X2 + Ti. (PW^) 
Jk. 

Here mi and 777-2 = 1 — mi are two constants that describe the mass distribution between the 
peaks, and as long as X2 is confined to the stable region x > we can expect that each solution 
to (PW) defines via 7771^ + rn2Sx2{t) an (approximate) solution to the original model (FP). 

In what follows we denote the width of the first peak by w{t) and aim to derive an asymptotic 
formula for the evolution of 'w{t) that involves only the dynamical constraint i{t). For simplicity 
we assume again that the data at time to are localized and well prepared in the sense of 

^>(x, to) ~ ^a;i(to)(^)' w{tQ) l^, -X»» < Xi(to) < -a;*, X* < ^2(^0) < a:**. 

A key ingredient to any asymptotic analysis of the widening phenomenon is to give an appropriate 
description of the position of the first peak. Our ansatz is to define xi{t) as solution of the 
characteristic ODE, that means we set 

TXi = a- H'{xi). (28) 
This ansatz has the following advantages. As long as the first peak is narrow, we have 

a « miH'ixi) + 7772-ff'(x2), ^ ~ 777lXl + m2X2, 

and hence xi and X2 evolve according to the two-peaks model (TP). Consequently, for < r ^ 1 
we can describe xi and X2 in terms of the quasi-stationary approximations of (TPM), whose 
dynamics is completely determined by the dynamical constraint £{t) and the mass distribution 
between the peaks, see (27). A further advantage of (28) is that it gives rise to a quite simple 
evolution law for w(t). 
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3.2.2 Formula for the width of the peak 

In order to analyze the growth of w{t) we introduce the rescaUng 

where both the spatial scaling factor X(t) and the rescaled time 9{t) will be identified below. We 
then have w{t) = iyX{t)W{9{t)), where W{e) is the width of R at e. 

We readily verify that (PWi) and (28) imply that the evolution of R is governed by 

redeR = + l,d^,R+^-^d,iyR). (30) 

For the subsequent considerations we now assume that Xi{t) is a given continuous function with 

3^1(^1) = "Xit and 

—X** < Xi{t) < —X* V to < t < ti , —X* < Xi(t) < X* V ti < t < ^3 

for some t^. Heuristically, is the time at which both peaks merge (continuously or discontinu- 
ously) according to the two-peaks model (TP). 

As long as w{t) is small - but possibly much larger than u - we can expand the nonlinearity 
according to 

H'{xi - u\y) - H'{xi) = H"{xi)y + 0{u\y) = H"{xi)y + 0{w). 
Neglecting the higher order terms and defining A and 6 as solutions to 

tX = -H"{xi)X, t9 = X{to) = 1, 9{to) = 0, 

the nonlinear PDE (30) transforms into the heat equation 

deR = d^R. 

Consequently, for large 9 the rescaled profile R evolves in an almost self-similar manner, that 
means we can approximate 

1 / y' 



Hfe,,)»_exp^-^J. (31, 

Notice that this approximation implies W{9) ~ for 9^1 and holds as long as W{to) is of 
order 1. In order to characterize the width of the original peak it remains to understand how A 
and 9 depend on t. A direct computation yields 

A(i) = exp (-^ ^* f{s) ds^ , 9{t) = i ^* exp (^^ j'^ f{s) ds^ di, 

where / abbreviates 

m ■.= H"{x,{t)). 

Due to < r < 1 and 

f{t)>0 VO<t<ti, 9{t)>0 V0<t<t3, 

we now infer that 9(t) ^ 1 for all to < t < t^. This implies w{t) = vX{t)-\j9{t), and hence 

w{tf = - Texp (-- f f[s)ds--\ dl (32) 



where a is the scaling parameter from (24). 

We have now identified an explicit formula for w(t), which involves only the function xi(t). 
Notice that for small r this function can be computed by the quasi-stationary two-peaks approx- 
imation, whose evolution is independent of r and completely determined by ra\ and l{t). 
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3.2.3 Asymptotic description of the widening 



Formula (32) can be restated as 



w{t) 




(33) 



where the function 



m := tf{s)ds 



(34) 



tl 



is non-positive on the interval (0, t^) since ti is the only root of / in this interval. We now aim 
to show that this formula implies that the critical time t2 is determined by the conditions 



More precisely, if (^(f) attains the value —a for some time t2 with < ^2 < ^3, then in the limit 
r ^ the width of the unstable peak explodes instantaneously at t2- If, however, <p{t) is larger 
than —a for all times ti < t < ts, then the width of the peak is exponentially small in r even at 
time ts. 

In order to prove these assertions we now derive rough estimates for w{t). At first we consider 
to < t < tl. In this case we have < ^{t) for all < i < t, so (33) gives 



Now let t > tl but suppose that t is sufficiently small such that (p{t) + a = 6 > 0. Using 0(i) < 
we then we estimate 



and find again that the width of the peak is exponentially small in r. Finally, we consider a time 
t with tl < t < ts and 0(i) + a = —6 < 0. To derive a lower bound for 'w{t) we now employ the 
continuity of / at ti as follows. For < e < t — ti and ti < i < ti + £ we estimate 



Hh) + a = 0, 



tl < t2. 



(35) 





<P{t}>f{ti){i-ti)-o{e) 



to find 




dt > 





Combining this with (33) gives 




and we conclude that the width of the first peak is exponentially large in l/r. 
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We finally emphasize that the equation for w(t) can be simplified for t ^ t2 as follows. Exploiting 
the continuity of / at t2 we find ^(t2 + ts) = —2a + f{t2)sT + o{sT). For each s we therefore have 



{t2 + sT) = {w\t2) + 0{u^))exp{- f{t2)s + o{l)), 



where o(l) means arbitrary small for small r. In particular, choosing i2 with |t2 — = o{l) such 
thati(;2(i2) =1 we find 



w 



'(t2 + sr) = (1 + 0(1)) exp (/3s), = -/(ts) > 0. (36) 



3.2.4 The mass splitting problem 

As explained above, at the critical time t2 ~ t2 we expect that the system undergoes a rapid 
transition from the unstable-stable configuration to a new stable-stable configuration. In order 
to describe this transition, in particular, to predict the mass distribution between the emerging 
stable peaks, we propose to study a simplified mass-splitting model (MS), which describes (PW) 
on the rescaled time scale s = (t — ^2) /r in the limit u ^ 0. In other words, (MS) consists of the 
equations 

dsg = d4{H'{x)-a)g), (MSi) 

X2 = a-H'{x2), (MS2) 



ds 



a = mi / H' {x)Qdx + m2X2, (MS3) 
Jr 



which have no diffusion and satisfy the constraint via -^{mixi + 1112X2) = 0. Notice that this 
equation is the Wasserstein-gradient flow for the energy 



£ = T-L = m\ I H (x) g dx + m2H {X2) ■ 



In view of the above discussion we now impose asymptotic initial conditions at s = —00, which 
reflect that the mass splitting process starts in an unstable-stable configuration and that the 
unstable peak is a rescaled Gaussian due to the entropic randomness. To this end we identify 
t2 = t2 and denote by xi and X2 the quasi-stationary two-peaks approximation at t2, i.e. we have 

xi = Xo{a), X2 = X+{a), mix i + 1712X2 = i{t2). 

In accordance with (29), (31) and (36) we now require that 

g{x,s) ^ ^ ^ .0 ^ exp P = -H"ixi)>0 (37) 

' ^ 2V7rexp(/3s) 4exp(2;3s)y' ^ ^ ^ ' 

weakly* in the sense of probability measures, and that 

X2is) X2. (38) 



The gradient flow structure of (MS) implies that each solution approaches a stable-stable config- 
uration in the limit s ^ 00. This reads 

X2{s) > X2, Q[x, s) > (mi - mi2)()xi +mi20x2> 
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where — x** < xi < — and x^.. < X2 denote the positions of the emerging stable peaks and 
< mi2 < mi is precisely the amount of mass that is transferred during the splitting process 
from the unstable region into the stable region x > x*. Of course, the asymptotic data at s = +00 
must comply with 

H'{xi) = H'{x2), {nil - iTii2)xi + {m2 + mi2)x2 = £ = niixi + m2X2, 

but these conditions do not determine the three quantities mi2, xi, and X2 completely. This is not 
surprising and reflects that the amounts of mass that are transferred towards the stable regions 
depends crucially on the asymptotic shape of the unstable peak. In our case, however, this shape 
is a rescaled Gaussian and therefore we expect that the data at s = +00 are uniquely determined 
by the data at s = —00. In particular, we conjecture that there is a unique and continuous mass 
transfer function M such that 

mi2 = M(mi, a). (39) 

Both the existence and continuity of M are not obvious because the mass splitting problem 
involves two subtle limits. At first one has to show that the asymptotic condition (37) gives rise 
to a well-posed initial value problem at s = —00. Second, one has to guarantee that solutions 
do not drift as s — )• 00 along the connected one-parameter family of equilibrium solutions. A 
rigorous justification of the mass splitting function M is beyond the scope of this paper and left 
for future research. 

It is, however, possible to compute numerical approximations of M using the method of 
characteristics. More precisely, restricting to N characteristics ^n{s) with n = 1...A'', the mass 
splitting problem (MS) can be approximated by 

^^„(s) = a{s) - H'iUs)), ^X2(s) = ais) - H'{x2{s)), (40) 



with 



N N 
"^1 rw//- / N\ , TTi, / , rni 



^(^) = +m2H'{x2{s)), i=^Y.^n{s) + m2X2{s) = const. (41) 

n=l n=l 

Moreover, to mimic the asymptotic initial conditions (37) and (38), we choose a small parameter 
£ and set 

en(0) = xi+ eeir^n/N), ^2(0) = X2, 

with crf^^ being the inverse of erf(^) = Tr~^^'^ exp {—x'^) dx. The resulting finite-dimensional 
initial value problem can be integrated numerically, for instance by means of the explicit Euler 
scheme, which satisfies the constraint i = const up to computational accuracy. In the limit 

. s s^oo . /- / \ s^oo ( xi for 1 < n < A''i2, 
X2(.) > X2, Us) > [ Ni2<n< N. 

The critical index A^i2 finally determines mi2 via mi2 = miNi2/N. 

The strategy for the numerical computation of M is therefore as follows. We sample the 
two-dimensional parameter space of (MS) and solve for each choice of the parameters the ODE 
(40) with fj as in (41) on a sufficiently large time interval. The results arc illustrated in Figure 
8 for the potential (2), where < mi < 1 and a = H'{xi) are, as above, regarded as the two 
independent parameters. 
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i^ii/i^i versus mi mulmi versus ct-/o-. 




Figure 8: Numerically computed values of mi2 for potential (2) and several values of toi and a. 



3.3 Limit dynamics: Switching, splitting, and merging events 

Combining the arguments from the previous sections we are now able to characterize the slow 
reaction limit of (FP) in terms of singlc-pcak and two-peaks configurations. To describe their 
evolution we consider piecewise constant mass functions m\{t) and m2{t) along with piecewise 
continuous position functions and X2{t), where we allow xi and X2 to be undefined on 

intervals with mi = and 1712 = 0, respectively. These functions are coupled by the constraints 

H' (xi) = H'{x2) = cr , I = niixi + 1712X2 , mi + 777-2 = 1 (42) 

with 

Xl < X2 , < mi, 7772 < 1 , 

and define the limit energy via 

E := m\H{x\) + m2H{x2)- 

To formulate the limit model wc first specialize to strictly increasing dynamical constraints and 
discuss possible generalizations afterwards. 

3.3.1 Limit model for strictly increasing constraints 

Sections 3.1 and 3.2 reveal that the slow reaction limit for constraints with (4) can be described as 
follows. The initial dynamics is governed by the quasi-stationary evolution of a single stable peak 
within the stable region x < — and hence it seems natural to define 7771 (t) = 1, which implies 
xi{t) = i{t) and renders X2 to be undefined. At the first switching time, which is defined by 
Xl = i = —X*, the single peaks becomes unstable and its width starts to widen. It may happen 
that the motion of the unstable peak is faster than the widening, and then the peak remains 
localized until it leaves the unstable region via i = Xj^. In this case we find a Typc-I phase 
transition as there is no splitting of unstable peaks; notice that this happens if £ is sufficiently 
large or if the scaling parameter a is sufficiently large. 

For Type-II transitions, however, the width of the peak becomes eventually large within 
the unstable interval. This means there exists a first splitting time with — x* < £ < x^ at 
which the system jumps from the unstable single-peak configuration to a stable-stable two-peaks 
configuration, where the masses and the positions of the emerging stable peaks is determined 
by the mass splitting problem. In particular, after the first splitting event the positions Xi, X2 
and the masses 7771, 7772 are well defined, and both stable peaks move according to the dynamical 
constraint until the first peak located at xi becomes unstable at the next switching time defined 
by Xl = Xit. Afterwards we must carefully investigate the unstable-stable evolution in order to 
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decide whether there is a further spUtting of the unstable peak or whether both peaks finally 
merge continuously or discontinuously. 

In summary, the limit dynamics can be described as illustrated in Figure 5, that means 
intervals of quasi-stationary transport are interrupted by singular times corresponding to the 
following types of events: 

Switching: The peak at xi enters the unstable region. 

Splitting: The unstable peak at xi splits and the system jumps to a new stable- 
stable configuration with decreased mass mi and increased mass m2. 

Merging: The peaks in an unstable-stable configuration merge either continuously 
(with xi = X* = X2) or discontinuously (with xi < x^ < X2), or there 
is only a single peak with mi = 1 that leaves the unstable region (with 

More precisely, each limit trajectory comprises K + 1 switching times, K splitting times, and a 
final merging time 

< tsw,0 < isp,l < *sw,2 < ■•• < tsp,K < tsw,K < ^me < ■ 

where we have K = and if > 1 for Type-I and Type-II transitions, respectively. Here < 00 

is defined by — •^**) so proper two-peaks configurations can only exist for t <t^. 

Notice that an infinite number of switching and splitting events is not possible because the 
splitting condition (35) implies a lower bound for isp,fc — isw,fc-i via 

a = / \H"{xi{t))\ dt < C{tsp,k - tsw,k) 

with C = sup|j.|<j.^ \H"{x)\. Notice also that we must truly book keep the switching events since 
(35) doe not determine tsp,fe but only tsp,fe — isw,fe-i- 

We now describe the flowchart from Figure 5 in greater detail. 

Intervals of transport. Between consecutive singular times, and likewise initially for < t < tsw,0) 
the dynamics of xi and X2 is governed by a rate-independent system of non-autonomous ODEs. 
More precisely, differentiating (42) with respect to time yields 

mi = m2 = 0, xi = — r^, X2 = r^, (43) 

Z(mi, 7712, Xl, X2) 7722, Xl, X2) 

where we recall that x\ and X2 remain undefined for m\ = and m2 = 0, respectively. These 
ODEs imply 

,,,, 

Z(mi, m2, Xl, X2) 

with 

Z{mi, m2, Xl, X2) = miH"{x2) +m2H"{xi). 

Since the dynamical constraint i{t) is given, the initial value problem to (43) can - at least in 
principle - be integrated. 
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Jump conditions at singular times. At switching times all functions are continuous, that means 
we have 

|[mi|(tsw,fe) = |[a:il(tsw,fc) = IHI(isw,fe) = |£^]|(tsw,fc) = 0, i = 1,2, 

where := f{t + Q) - f{t-0) with f{t±0) = \ims\of{t±d) denotes the jump of the 

function / at time t. 

At a splitting time all quantities do jump, where the jump heights are determined by the mass 
splitting problem (MS). In particular, the gradient flow structure of (MS) ensures that 

m{tsp,k)<o, 

and using the mass transfer function M from (39) we find 

|[?7i2]|(tsp,fe) = -\[mi]\{tsp^k) = M(mi(tsp,fc - 0), o-(tsp,fe - 0)). 

The jumps of xi, X2 and a are then determined by the algebraic constraints (42), and can be 
used to reinitialize the ODEs (43) and (44). 
Finally, at the merging time we have 

Xl(tme + 0) = X2(tme + 0) = ^(tme) , (^(Ue + 0) = i7'(^(tme)) , 

and the concavity condition (A3), see the discussion at the end of Section 3.1.2, provides 

|[^l(tmc) <0 

with strict inequality for discontinuous merging. After the merging the precise values masses mi 
and m2 are actually undefined, but it seems natural to set 

TOi(tme + 0) = 0, m2(tme + 0) = l. 



Determining the next singular time. After the k^^ splitting time (we set tsp,o := to describe the 
initial evolution), the system is in a stable-stable configuration and the subsequent switching time 
^sw.fc is the smallest time larger than 4p,fe such that cr(tsyf^k) = (^*- This implies xi(tsw,fe) = —x*, 
X2{tsw,k) ~ 3^**, and 

^(*sw,fe) = -mix* + m2X**, 

which determines tsw,fc uniquely since mi and m2 are known. 

After the switching event at tsw,fe-i we have to decide whether the next singular time corre- 
sponds to splitting or merging according to the conditions 

discontinuous merging: Z{mi, m2, xi(tme), a^2(imc)) = 0, 
continuous merging: xi{t 



splitting: / H" [xi {t)) dt + a = 0. 

*sw,fc-l 



To discuss this decision in a simple case, let us suppose that i is constant in time. Due to d < 
we can replace t by a, and direct computations provide the following nonlinear equations 

discontinuous merging: Z{mi, m2, ^o(c"me)5 ^+{'^me)) = 0, 
continuous merging: a^e = — cr*. 



J a. 



Zjmi, m2, Xo{a), X +{a)) 



splitting: I ' t^..,^ /~x\ do" = at. 
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In particular, both merging conditions are rate independent, whereas the switching condition 
depends on i and cannot be satisfied for large I. This is not surprising because the merging 
conditions are completely determined by the quasi-stationary two-peak approximation, whereas 
splitting only happens if the peak widening due to the separation of characteristics is faster than 
the transport due to the dynamical constraint. 

Instead of solving (43), we can now vary a, starting from cr* and moving towards —a*, and 
check which conditions is satisfied at first. In any case, the time tev of the next singular event 
can be computed by 

, _, _ r* Zimi,m2,Xo{a),X+{a)) 

iev H"iXoia))H"{X^{a)) 

3.3.2 Limit model for non-monotone constraints 

We finally derive an alternative description of the slow reaction limit that covers arbitrary dy- 
namical constraints i. Decreasing constraints can produce stable-unstable configurations with 
xi < —X* and —x^ < X2 < x^t, which in turn can split or merge. However, these effects can 
easily be described by adopting the formulas from sections 3.1 and 3.2. The truly new effect is 
that non-monotone constraints can trigger inverse switching events, that means, for instance, an 
unstable-stable configuration can become stable-stable via xi = — a;* with xi < 0. 

For general constraints, it is convenient to describe the slow-reaction limit in terms of the 
following variables: (i) the dynamical multiplier cr, (m) an internal variable <p G [—a, 0] as intro- 
duced in Section 3.2 to control the width of an unstable peak, and (Hi) three nonnegative masses 
m_, mo, and m+ with 

^ mi = 1 , JJ mj = 

ie{-,o,+} ie{-,o,+} 

to describe the several types of single-peak and two-peaks configurations. Using these parameters, 
the peak positions are dependent variables which satisfy Xi = Xi{x) for i G {— ,0, -|-} as long as 
Xi is well-defined. 

configuration masses rrii range for a range for (f) further constraint 



5_ 


m_ 


= 1 


— oo < a < -l-CT* 


<j> = 









s+ 


m+ 


= 1 


— cr* < a < +00 


<!> = 









So 


mo 


= 1 


— cr* < a < -|-(T* 


—a < 4 


)<0 






T-+ 


mo 


= 


— cr* < a < 


cf> = 









T-o 


m+ 


= 


— (T* < 0" < -|-(T* 


—a<4 


)<0 




> 


%+ 


m_ 


= 


—a<^ < a < -\-a^: 


-a < 4 


) < 


Zo+ 


> 



Table 2: List of all possible single-peak configurations <Sj and two-peaks configurations Tij. 

At each non-singular time, the system is confined to one of the sets defined in Table 2, where 
Z_o and Zo+ are abbreviations for 

Z_o := m_^o(cT) + moA_{a) , Zo+ := mo^+(o-) + m+Ao{a) , 

and the functions Ai are defined by 

Ai{a):=H"{Xi{a)), ie{-,0,+}. 
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Within each of these sets, the peaks are transported according to the quasi-stationary two-peaks 
approximation and the widening of unstable peaks is governed by (f) as described in (33) and (34). 
This reads 

rhi = 0, a = 1^ miAi{a), ^ = X{mo^o}^o(o-) • (45) 

ie{-,o,+} 

On the other hand, singular events happen when the system reaches the boundary of either Si 
singular event condition possible jumps 



switching 


a = 


+0"* 


S- ^ So, 


T-+ ^ To+ 




a = 


— fj* 


S+ ^ So, 


T-+ ^ T-o 


inverse switching 


a = 




So ^ S-, 


To+ ^ T-+ 




a = 


— cr* 


So ^ S+, 


T-o ^ T-+ 


splitting 


4> = 


—a 


So 71+, 


T-o ^ T- 


discontinuous merging 




= 


T-o ^ 5_ 






Zo+ 


= 


%+ 1-^ s+ 




continuous merging 


a = 


+cr* 


T-o ^ 5- 






a = 


— cr* 


To+ 1-^ S+ 





Table 3: List of all possible singular events. 



singular event subcase behavior of variables 



switching a = m_ ■<->■ mo 

a = — cr* mo -H- m+ 

inverse switching a = m_ -H- mo, (f> 

a = — cr* mo ■<->■ m+, ^ — >■ 

splitting mo 0, |[m_]| > 0, |[m+]| > 0, (f) ^ 0, \[a]\ = 

merging into 5+ m+ — )■ mo -|- m+, mo — )■ 0, — > 0, < 

into S- m_ — >■ m_ -|- mo, mo ^0, 0^0, |[(t]| > 

Table 4: Update rules for variables m,, ct, and at singular events. 

or Tj. More precisely, examing the conditions for switching, inverse switching, splitting, and 
merging we arrive at list from Table 3, where we now interpret trivial continuous merging as 
inverse switching of single-peak configurations. The corresponding jump and update rules for the 
variables mi, a, and are summarized in Table 4. 

The slow reaction dynamics of two-peaks initial data can now be integrated iteratively by (i) 
moving via (45) along either one of the sets Si and Tij, and (ii) jumping to another set when 
reaching the boundary. Notice that at both splitting and discontinuous merging events the system 
jumps to inner points and that splitting events require to solve the mass splitting problem. 
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